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1  INTRODUCTION 


Although  the  present  rate  of  increaise  of  CO2  in  the  atmosphere  is  well 
documented,  changes  in  the  next  century  cannot  be  accurately  predicted 
because  the  CO2  budget  is  so  poorly  understood.  The  atmospheric  con¬ 
centration  is  increasing  1.4  ppm/year,  corresponding  to  an  increase  in  the 
atmospheric  inventory  of  3  Gt/year  {I  Gt  is  10^  metric  tons  or  10'^  kg). 
This  observed  increase  is  less  than  half  of  the  estimated  8.7  Gt/year  released 
by  all  sources  to  the  atmosphere  [10].  Most  of  the  residual  has  been  assumed 
to  enter  the  ocean  [15].  This  cannot  be  verified  directly,  however,  since  the 
total  dissolved  inorganic  CO2  in  the  ocean  is  so  much  larger  than  the  cu¬ 
mulative  release  of  fossil  fuel  since  1850  that  the  change  due  to  fossil  fuel 
is  small  compared  to  natural  variability.  As  a  result,  global  CO2  budgets 
are  attempted  by  calculating  the  fluxes  between  ocean  and  atmosphere  and 
between  land  and  atmosphere. 

There  are  two  approaches  to  CO2  fluxes  between  ocean  and  atmosphere. 
The  first,  widely  used  in  geochemical  studies,  uses  bulk  formulas  based  on 
the  difference  of  the  partial  pressure  of  CO2  denoted  PCO2,  between  air  and 
water,  and  on  Uio,  the  wind  speed  10  m  above  the  sea  surface.  The  bulk  ap¬ 
proach  is  not  well  established,  however,  and  differences  exist  in  the  formulas, 
in  methods  of  measuring  partial  pressures,  and  in  the  spatial  and  temporal 
coverage  of  the  data  used  to  evaluate  the  formulas.  Consequently,  global 
calculations  of  net  CO2  fluxes  also  vary  widely.  The  second  approach  to  CO^ 
fluxes,  developed  by  boundary  layer  turbulence  investigators,  is  more  recent 
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and  utilizes  new  fast-responding  detectors  for  CO2  concentration.  These  are 
necessarily  point  measurements,  but  similar  direct  measurements  of  heat  and 
water  vapor  fluxes  have  been  used  to  develop  and  validate  bulk  formulas  for 
those  variables.  Although  the  instrumentation  is  still  being  developed,  the 
results  seem  consistent  with  the  more  traditional  measurements  of  water  va¬ 
por  and  heat  fluxes.  The  direct  COj  measurements,  however,  give  fluxes 
much  larger  than  bulk  formulas  used  by  the  geochemists,  producing  vigorous 
controversy. 

Before  considering  how  the  flux^  are  estimated,  the  treatment  by  Broecker 
and  Peng  [14]  is  summarized  to  provide  some  background  about  the  chem¬ 
istry  of  carbon  dissolved  in  the  sea. 
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2  CARBON  CHEMISTRY 


C&rbon  dioxide  gas  crosses  the  sea  surface  primarily  because  of  the 
reaction 


HjO  +  COj  +C0j  ^2HC0J,  (2-1) 

where  CO3  is  the  carbonate  ion  and  HCOj  is  the  bicarbonate  ion.  Without 
this  reaction,  CO2  would  exist  in  the  upper  ocean  only  ais  dissolved  gas,  and 
the  ocean  would  take  up  only  3%  of  the  fossil  fuel  CO2  [2].  The  concentrations 
«ire  related  by 


K'  = 


[HCO3-P 


(2-2) 


(Ha0][C0^](CO2]’ 

where  the  equilibrium  constant,  K\  is  a  function  of  salinity,  temperature,  and 
pressure.  Dissolved  inorganic  carbon  occurs  in  the  upper  ocean  in  the  three 
reduction  states  of  CO2,  and  less  than  1%  of  the  total  dissolved  inorganic 
carbon,  denoted  ECO2,  is  found  in  dissolved  CO2  gas.  Consequently,  the 
concentration  of  total  dissolved  inorganic  carbon  can  be  approximated  as 
the  sum  of  its  two  major  components 


[ECO2]  =  [HCOJ]  +  [COj] 


moles 

.  . 


(2-3) 


Owing  to  the  buffering  mechanism  expressed  by  Equation  (2-1),  the 
time  for  CO2  in  surface  water  to  equilibrate  with  the  atmosphere  is  about 
10  times  longer  than  for  other  gases — to  change  PCO2,  it  is  also  necessary  to 
change  the  mn.ch  larger  inventories  of  CO3  and  HCO3 .  Broecker  and  Peng 
[2]  estimate  the  step  response  time  of  pC02  and  SCO2  in  surface  water  as 
about  one  year. 
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Reflecting  the  important  role  of  the  negative  ions  COj  and  HCOj  in 
the  balance  of  electrical  charge  in  the  sea,  alkalinity  is  defined  as  the  excess 
of  positive  charges  to  be  balanced  by  these  ions 


[Alkj  =  [HCO3-I  +  2[C0rl 


moles 

.  kg  . 


(2-4) 


Therefore,  alkalinity  is  also  the  sum  of  the  concentrations  of  positive  charges 
less  the  concentrations  of  the  remaining  negative  ions. 


Because  the  activity  of  water  is  identically  one,  Equations  (2-2),  (2-3), 
and  (2-4)  can  be  manipulated  to  give  the  concentration  of  carbon  dioxide  gas 
in  terms  of  the  concentration  of  SCO2,  adkalinity,  and  K' 


[CO2]  = 


I  (2(SC02]  -  [Alk])^ 
K'  [Alk)-[EC02] 


moles 


I  kg  J 


(2-5) 


Consequently,  one  method  of  determining  (CO2]  is  to  measure  (ECO2]  and 
[Aik],  and  use  the  known  dependence  of  K'  on  salinity,  temperature,  and 
pressure.  (The  concentration  can  be  converted  to  mctss  units  using  44.01  g 
for  the  molecular  weight  of  CO2.  For  example,  a  concentration  of  1  /xmole 
CO2  per  kilogram  of  seawater  corresponds  to  0.044  mg  CO2  per  kilogram  of 
seawater.) 


In  air,  the  partial  pressure  of  CO2  is  the  total  pressure  of  the  air  times 
the  mole  fraction  of  CO2  in  the  gas  mixture.  In  water,  pCOa  is  the  partial 
pressure  of  CO2  in  air  that  is  in  equilibrium  with  the  solution  at  the  in-situ 
temperature  and  pressure.  Thus,  the  CO2  concentration,  [C02],  is  related  to 
the  partial  pressure,  PCO2,  by 

PCO2  =  [atm],  (2  -  6) 

Q 
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where  o  is  the  solubility  of  CO3  in  water,  with  units  of  moles  COj  per 


kilogram  atmosphere.  Combining  with  Equation  (2-5), 


1  (21£C0,1  -  lAlkp’ 
"  aK'  (AUtl  -  (ECO,) 


[atm|. 


(2-7) 


As  seen  in  Figure  1,  over  the  range  of  sea  surface  temperatures,  a  decreases 
by  about  a  factor  of  two  with  increasing  temperature.  (Weiss  [23]  gives  a  in 
terms  of  fugacity;  Heimann  and  Monfray  [8]  express  Weiss’s  results  in  terms 
of  partial  pressure.)  Owing  to  the  temperature  dependence  of  both  q  and 
K',  their  product  varies  by  a  factor  of  3.7  between  0°C  and  30°C  (Figure  2). 
As  a  result,  the  equilibrium  partial  pressure  increases  with  increasing  tem¬ 
perature  by  4%/°C  when  total  carbon  and  alkalinity  are  constant.  Because 
the  difference  in  partial  pressures  between  air  and  water  is  often  only  a  few 
percent,  the  temperature  dependence  of  aK*  greatly  complicates  estimates 
of  the  global  rate  of  CO2  exchange. 

Ray  Weiss,  at  Scripps  Institution  of  Oceanography,  measures  ApC02 
by  drawing  samples  of  air  and  water  while  a  ship  is  underway.  The  partial 
pressure  of  the  ocean,  (pC02)oce?  is  determined  by  equilibrating  the  water 
and  air  samples  at  the  in-situ  temperature  and  pressure.  The  samples  are 
then  dried  and  analyzed  in  a  gas  chromatograph,  and  the  results  are  reported 
in  terms  of  moles  of  CO2  per  mole  of  dry  air.  Examples  are  shown 

in  Figures  3  and  4.  Conversion  of  requires  knowledge  of  the  humidity 
expressed  eis  the  partial  pressure  of  water,  phjO,  in  terms  of  which 


pCOi  =  (/\ot«i  -  PHjo)a:c02 


(2-8) 


where  Ptotai  is  the  total  pressure  (this  is  equation  10  of  Weiss  and  Price  [25]). 
To  determine  the  average  flux,  it  would  be  appropriate  to  use  the  average 


0  2.000  4,000  6,000  8,000  10,000 


Distance  (km) 


Figure  3.  Mole  fractions  of  CO^  in  air  and  water  measured  continously  during  Leg  12 

of  the  INDOMED  cruise  from  Puerto  Rico  to  Buenos  Aires  from  2 1  September 
to  1  November  1978  [Reference  24].  As  the  differences  in  mole  fractions  are  dose 
to  differences  in  partial  pressures,  the  large  excess  of  ACO^  in  the  water  near  the 
equator  reflects  strong  outgassing  up  the  upwclling  water. 
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TTO/NAS  Leg  b 


Distance  (km) 


Figure  4.  Mole  fractions  of  CO^  in  air  and  water  measured  continously  during  Leg  5  of  the 
North  Atlantic  Study  of  the  Transient  Tracers  cruise  from  Glasgow  to  Rcyjavik 
between  19  June  and  15  July  1981  [Reference  24],  The  large  deficit  of  AGO,  reflects 
strong  absorption  in  the  North  Atlantic  at  high  latitudes. 


9 


(PHjo)air)  since  xcoa  ‘s  relatively  constant  over  time.  As  evident  in  Figure  5, 
the  correction  for  PUjO  is  small,  but  it  greatly  increases  with  temperature. 

Biology  is  also  an  important  aspect  of  the  carbon  chemistry  of  the  sea 
and  one  that  is  likely  to  be  strongly  correlated  with  wind  speed.  Deepening 
surface  mixed  layers  often  entrain  nutrient-rich  waters  from  the  thermocline, 
producing  an  explosive  growth  of  phytoplankton  in  the  mixed  layer.  These 
issues  are  beyond  the  scope  of  this  preliminary  survey,  but  are  important  and 
need  careful  study. 


10 


T  - - -  ' - '  ' - ^ 

5  10  15  20  25  30 

°C 


Figure  5.  The  partial  pressure  of  water  vapor  in  air  is  a  small  correction,  but  increases  strongly 
with  temperature  [Reference  25]. 
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3  BULK  FLUXES 


The  bulk  flux  is  formulated  as 

FcOi  =  E  •  ((pC02)*tin  -  (pC02)oce]  (3-1) 

=  E  •  ApCOi  [mol  CO2  s"*],  (3-2) 

where  E  is  the  gas  exchange  coefficient  with  units  of  moles  CO2  m“^  s~*  atm"’ , 
and  ApCOi  is  in  atm. 

The  atmosphere  is  well  mixed  in  CO2,  with  a  mean  concentration  close 
to  350  ppm.  The  principal  variability  is  an  annual  cycle  of  about  5  ppm  (re¬ 
sulting  from  the  yearly  change  in  plant  respiration)  and  a  mean  north/south 
gradient  of  about  3  ppm  (resulting  from  greater  release  of  fossil  fuel  in  the 
northern  hemisphere).  In  surface  waters,  oceanic  concentrations  range  from 
150  to  450  ppm.  Therefore,  changes  in  the  driving  potential  for  CO2  exchange 
acruss  the  sea  surface  result  largely  from  changes  of  (pC02)oce- 

The  partial  pressure  of  CO2  in  the  surface  water  is  measured  by  two 
methods.  One,  described  earlier,  is  Weiss’s  which  uses  an  air/water  equilibra- 
tor  connected  to  a  gas  chromotograph.  This  approach  operates  continuously 
while  the  ship  is  underway  and  determines  (pC02)4tm  and  (pC02)oce-  The 
other  method,  used  by  Taro  Takahashi  at  Lamont,  is  to  collect  bulk  water 
samples,  determine  the  total  carbon  and  pC02  at  a  standard  temperature, 
and  then  correct  pC02  to  in-situ  values  using  the  temperature  dependences 
of  the  solubility  and  the  equilibrium  constants.  Comparisons  of  simultaneous 
measurements  reveal  root  mean  square  (rms)  differences  of  6-7  //atm  between 
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the  two  techniques,  and  along  some  ship  tracks  biases  of  2-3  patm  have  been 
found  [24].  Although  these  errors  are  relatively  small  compared  to  the  con¬ 
centration  of  CO2,  they  sometimes  exceed  the  difference  in  concentration 
between  air  and  water  over  large  areas  of  the  ocean.  For  the  first  4,000  km  of 
Weiss’s  measurements  from  San  Juan  to  Buenos  Aires  (Figure  3),  the  differ¬ 
ence  in  concentrations  of  COj  between  air  and  water  is  no  more  than  10  ppm 
and  changes  sign  several  times.  In  other  locations  the  errors  will  indeed  be 
small  relative  to  the  contrast  in  partial  pressures.  For  instance,  upwelling 
produces  a  large  excess  of  CO2  in  the  surface  water  at  the  equator  causing 
outgassing.  Another  caise  of  large  contrasts  occurs  at  high  latitudes  during 
summer  in  the  North  Altantic  (Figure  4).  For  several  thousand  kilometers, 
the  oceanic  concentration  is  about  130  ppm  less  than  the  atmospheric  value 
causing  strong  absorption. 

The  gas  exchange  coefficient  is  usually  expressed  as 

E  =  a  Vp  [moles  CO2  m~^  s~*],  (3  —  3) 

with  Vp  as  the  gas-exhange  velocity.  Also  known  as  the  piston  velocity,  Vp  is 
not  well  known.  Presently,  the  most  realistic  parameterization  [11]  expresses 
it  in  terms  of  temperature  and  wind  speed 

'  (^)''''4.7xl0-^f/.o  t/,o<3.6 

Vp  =  -  (7.92  X  IQ-^Uio  -  2.68  x  lO*^)  3.6  <  (7,0  <  13 

(1.64  X  10-*f/,o  -  1.369  X  lO--*)  13  <  (7,0 

(3-4) 

where  Vp  is  in  m/s,  Uio  is  the  wind  speed  (in  m/s)  at  10  m  above  the  sea 
surface,  and  Sc  is  the  Schmitt  number  for  CO2  in  water. 


The  temperature  dependence  occurs  via  the  Schmitt  number,  defined 


as  the  ratio  of  the  kinematic  viscosity,  to  the  diffusivity  of  CO2  in  water, 
AccOj’  which  vary  with  temperature.  As  i/  a:  1  x  10"®  m^  s"‘  and 

=  1  to  2  X  10"®  m^  s"*,  Sc  as  600  at  20°C.  Between  0®C  and  30®C, 
5c/600  varies  more  than  a  factor  of  three,  decreasing  with  increasing  tem¬ 
perature.  (5c  also  varies  with  salinity,  but  the  effect  is  slight  and  can  be 
ignored  for  CO2  studies.)  Owing  to  this  dependence  of  Sc,  Vp  increases  with 
temperature  by  about  a  factor  of  two  (Figure  6). 

This  increase  of  Vp,  however,  is  offset  by  the  decrease  of  the  solubility,  a. 
As  a  result,  the  CO2  exchange  coefficient,  E,  varies  little  with  temperature 
(Figure  7). 

Liss  and  Merlivat  [11]  proposed  the  wind  speed  dependence  in  Equation 
(3-4)  based  on  an  approximate  fit  to  measurements  in  wind  tunnels  and  over 
lakes.  The  increase  of  V'p  with  wind  speed  (Figure  8)  differs  from  the  transfer 
coefficients  for  heat,  water  vapor,  and  momentum,  which  show  little  change 
with  Uio  when  the  atmosphere  is  weakly  stratified.  Transfer  coefficients  for 
other  quantities  increase  nearly  linearly  with  Uio,  but  gases  seem  to  have  a 
stronger  dependence.  This  is  not  surprising  since  the  molecular  diffusivity 
KCO2  small,  resulting  in  a  resistance  to  transport  almost  entirely  from 

the  water  phase. 

Bubbles  may  play  an  important  role  in  gas  transfer  for  two  reasons. 
First,  bubbles  greatly  increase  the  total  area  of  the  air-sea  interface.  Be¬ 
cause  the  resistance  to  transport  is  so  completely  in  the  water,  bubbles  lasting 
longer  than  turbulent  eddies  in  the  air  can  still  be  important  to  gas  transfer. 
Second,  the  global  average  of  ApC02  is  of  the  order  of  10  /ratm,  much  smaller 


Figure  6.  Temperature  dependence  of  the  piston  veloaty,  where  Sc  =  10’ 
[Reference  8]  based  on  data  in  [Reference  7]. 
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Figure  7.  Total  temperature  dependence  of  the  exchange  velodty.  £,  varies  little  with 

temperature,  owing  to  compensating  effects  in  the  solubility,  alpha,  and  in  the 
Schmitt  number  dependence  of  the  piston  velodty. 
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Figure  8.  Variation  of  the  piston  vclodtY  V .  with  wind  speed  (solid)  compared  with 
observations  and  other  models  [Reference  11]. 


than  (pC02)*tmj  which  is  approximately  350  /iatm.  Small  bubbles  are  rou¬ 
tinely  observed  several  meters  below  the  surface  [5]  and  their  concentration 
increases  rapidly  with  Uio,  perhaps  as  U^q.  Owing  to  surface  tension,  the 
pressure  inside  small  bubbles  is  higher  than  the  hydrostatic  pressure.  Thus, 
the  effect  of  bubbles  on  gas  transfer,  discussed  by  Thorpe  [22],  could  cause 
E  to  vary  rapidly  with  wind  speed.  In  a  similar  way,  small  ripples  which 
modulate  the  thin  diffusive  boundary  layer  below  the  air-sea  interface  may 
be  important  for  gas  transfer.  Because  much  of  the  information  about  the 
piston  velocity  comes  from  laboratory  wind  tunnels,  it  is  not  clear  how  well 
these  measurements  reflect  gas  transfer  in  a  fully  developed  wind-driven  sea. 

Owing  to  the  non-linear  dependence  of  E  on  t/jo,  instantaneous  wind 
speeds  are  essential  for  evaluating  Equation  (3-4).  Laboratory  and  local 
oceanic  field  studies  of  gas  transfer  are  based  on  instantaneous  winds,  and 
the  non-linear  dependence  does  not  permit  the  use  of  temporally  or  spatially 
averaged  winds.  As  a  minimum,  reasonably  accurate  probability  density 
functions  should  be  used  to  calculate  the  average  exchange  coefficient,  <  E  >. 
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4  TURBULENT  FLUXES 


Caxbon  dioxide  fluxes  axe  obtained  directly  in  the  atmospheric  boundary 
layer  by  measuring  the  correlation  between  fluctuations  in  vertical  velocity 
and  fluctuations  in  the  CO2  concentration 

®  ^  s“*],  (4-1) 

where  to'  is  in  m/s,  and  d  is  in  mg  C02/m^.  Turbulence  investigators  ex¬ 
press  Fqq^  in  units  of  mg  s,  in  which  case  a  =  1.  The  moleculcir 

weight  of  CO2  is  used  to  convert  to  geochemical  units,  i.e.,  a  =  {1/44  x  10^ 
mg  per  mole  CO2).  Some  geochemists  use  units  of  mole  CO2  m“^  year"’ 
[3].  These  can  be  converted  to  the  units  of  the  turbulence  investigators  us¬ 
ing  1  mole  CO2  year"’  =  1.40  x  10"^  mg  CO2  s"’.  Similarly,  CO2 
concentrations  can  be  converted  by  1  ppm  CO2  =  1-83  mg  CO2  m"^. 

Beginning  over  a  decade  ago,  turbulence  measurements  of  CO2  have 
developed  slowly  [9],  [12],  [17],  [26]  and  generated  sharp  controversy.  To 
the  extent  that  they  have  been  compared,  the  turbulent  results  exceed  the 
bulk  fluxes  by  factors  of  5  to  10.  Rejecting  the  turbulence  fluxes,  most  of 
the  geochemists  attacked  [3],  pointing  out  the  low  signal-to-noise  levels  of 
the  CO2  measurements  and  citing  agreement  between  their  fluxes  and  bulk 
constraints  imposed  by:  1)  distributions  of  naturally  occurring  between 
ocean  and  atmosphere,  2)  evasion  rates  of  radon  from  surface  waters,  and 
3)  global  distributions  of  bomb  ‘^C.  Smith  and  Jones  [18]  and  Wesley  [27] 
defended  their  measurements. 
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Caxbon  dioxide  fluxes  are  very  low  over  the  ocean  compared  to  those  over 
letnd,  where  growing  crops  produce  fluxes  of  1  to  2  mg  C02/m^  s  (2-4  x  10~® 
moles  COa/m^  s),  values  which  match  the  fluxes  predicted  by  agricultural 
measurements.  By  contrast,  0.083  mg  COj/m^  s  wais  the  largest  magnitude 
observed  by  Smith  et  al.  [19].  With  such  weak  fluxes,  noise  and  drift  are  key 
issues  for  evaluating  fluxes  over  water. 

The  turbulence  measurements  are  made  with  two  sensors  located  within 
approximately  1  m  of  each  other.  Sonic  anemometers  detect  w'  and  infrared 
absorption  detectors  sense  C02.  The  new  generation  of  sensors  detects  the 
absorption  from  a  strong  infrared  source  (at  a  temperature  of  about  1,000  K) 
over  a  path  in  the  air  of  0.2  to  0.25  m.  Attenuation  is  compaied  between 
a  wavelength  of  4.3  fim,  center  of  the  strong  CO2  absorption  bands,  and 
3.9  /im,  which  is  not  absorbed.  To  remove  the  effects  of  absorption  by  water 
vapor,  attenuation  is  also  compared  between  2.6  pm,  center  of  the  water 
vapor  absorption  bands,  and  3.9  pm.  Smith  et  al.  report  peak-to-peak  noise 
of  0.3  mg  C02/m^,  a  factor  of  ten  lower  than  for  their  previous  measurements. 

The  new  sensors  reveal  the  same  discrepancy  with  the  bulk  formulas  as 
before.  For  Uw  =  3  to  9  m/s  Smith  et  al.  [19]  report 

Kp  =  6.5  X  10-*  Uio  [m/s).  (4  -  2) 

Comparison  with  the  bulk  formulas  shows  that  the  direct  measurements 
are  larger  by  factors  of  5  to  10  (Figure  9).  Furthermore,  cospectra  of  w'  and 
c*  are  very  similar  to  those  of  w'  and  T'  and  of  w'  and  q\  where  q'  is  the 
fluctuation  in  water  vapor  (Figure  10). 
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Figure  10.  Cospectra  of  v/  with  (1)  CO^  concentration  (solid  circles),  (2)  temperature  (open 
circles),  and  (3)  water  vapor  {+)  [Reference  13].  The  structure  of  the  cospectrum 
with  CO  is  very  similar  to  that  of  the  other  variables. 


We  conclude  that  the  turbulent  measurements  seem  as  carefully  done 
as  those  for  other  atmospheric  variables  and  cannot  be  dismissed.  Neither, 


however,  can  we  dismiss  the  bulk  constraints  cited  by  the  geochemists  as 
justification  that  the  bulk  formulas  give  the  correct  magnitude  for  average 
fluxes  over  large  areas.  The  discrepancy  probably  results  from  differences  in 
the  space  and  time  scales  of  the  two  approaches,  as  argued  by  Smith  and 
Jones  [18].  Nevertheless,  the  controversy  needs  to  be  resolved  soon. 
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5  MEASURING  WINDS 


It  is  tempting  to  compute  time  and  spaee  averages  of  COj  exchange 
from  the  average  wind  speed,  and  this  is  frequently  done.  For  example,  Tans 
et  al.  [21]  use  a  linear  dependence  on  wind  speed 

0  Uio<0 

0.016(f/io  -  3)  3<t/,o 

which  they  attribute  to  Peng  and  Takahashi  [14].  To  compensate  for  the 
lineaLT  dependence,  this  expression  doubles  Vp  relative  to  the  Liss  and  Mer- 
livat  formulation.  If  Vp  and  E  varied  linearly  with  wind  speed  {as  they  do 
for  heat  transfer),  this  would  be  acceptable  so  long  as  ApC02  fluctuations 
were  constant  or  varied  with  wind  variations.  But,  as  the  empirical  law  in 
Equation  (3-4)  indicates,  the  relation  between  transfer  rate  and  wind  speed 
is  nonlinear.  It  is,  therefore,  unacceptable  to  replace  the  aA'erage  of  E  with 
E  evaluated  at  the  average  wind  speed,  i.e.,  to  replace  <  E{Uio)  >  with 
E{<  U\Q  >).  As  a  minimum,  a  reasonably  accurate  probability  density  of 
wind  speed  is  required  to  calculate  an  average  exchange  coefficient.  Beyond 
this,  the  covariation  of  ApCOj  and  wind  speed  could  be  deduced  and  used 
for  better  estimates  of  average  fluxes. 

Direct  observations  of  wind  speed  from  oceanographic  buoys  are  suf¬ 
ficiently  accurate  for  estimating  fluxes,  but  buoys  are  too  scarce  to  be  of 
much  help  for  global  distributions.  Routine  ship  observations  are  much  more 
numerous.  Many  are  based  on  visual  identification  of  the  sea  state,  which 
represents  an  average  of  an  hour  or  so,  owing  to  the  finite  growth  rate  of  the 


moles  CO2 


m^  year  patm 


(5-1) 
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short  waves  most  affecting  the  appearance  of  the  sea.  Because  the  observers 
are  frequently  well  trained,  these  data  are  consistent  and  useful  [16],  so  much 
so  that  scatter  between  such  wind  estimates  is  about  25%  and  averages  from 
different  observers  agree  to  about  10%  [20].  Well-placed  ship  anemometers 
with  good  exposure  appear  to  agree  with  buoy  observations  to  5%,  but  the 
same  ship  on  a  different  course  relative  to  the  wind  may  yield  errors  of  30%. 
The  main  problem  with  direct  observations,  however,  is  their  low  density, 
particularly  in  the  southern  hemisphere  and  in  high  latitudes  during  winter. 
Efforts  are  underway  to  improve  the  ability  of  operational  weather  forecast¬ 
ing  models  to  estimate  the  wind  information  needed  for  flux  estimates,  but 
success  has  not  been  demonstrated  where  data  are  sparse. 

Can  adequate  winds  be  determined  from  satellites?  Esse.itially  three 
types  of  satellite  sensors  can  be  used  to  estimate  wind  speeds.  Radar  scat- 
terometers  are  designed  to  measure  wind  speed  and  direction.  Comparisons 
of  the  Seasat  scatterometer  with  direct  wind  observations  from  oceanographic 
buoys  indicate  a  scatter  of  about  ±2  m/s  for  speeds  of  0  to  12  m/s.  Using 
the  same  phenomenon,  i.e.,  increasing  surface  roughness  with  wind  speed, 
radar  altimeters  measure  wind  speed  (but  not  direction).  Because  of  the 
vertical  angle  of  incidence,  these  measurements  should  be  slightly  less  ac¬ 
curate  than  those  from  a  scatterometer,  but  Chelton  and  McCabe  f4|  find 
the  degradation  to  be  slight.  Finally,  passive  microwave  sensors  also  respond 
to  wind  speed  variations,  presumably  because  wind-generated  surface  rough¬ 
ness  affects  the  microwave  emissivity.  Comparison  of  winds  derived  from  the 
SMMR  radiometer  on  Seasat  with  winds  from  buoys  showed  a  scatter  as  little 
as  ±1  m/s  [Ij. 
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Both  direct  and  satellite  wind  observations  have  instrumental  errors 
and  potential  biases.  Nevertheless,  by  comparing  Seasat  winds  with  direct 
observations  and  analyzed  field  reports,  Cheiton  [private  communication] 
finds  that  the  largest  errors  in  wind  statistics  come  from  inadequate  sampling. 
Even  though  a  satellite  makes  many  more  observations  than  all  the  ships 
at  sea,  the  influence  of  unobserved  wind  events  on  the  climatology  is  still 
the  greatest  source  of  error.  Thus,  the  best  estimates  of  wind  climatologies 
would  be  obtained  by  combining  all  satellite  sensors  with  all  available  direct 
observations. 

No  scatterometers  are  presently  operating,  but  the  European  satellite 
ERS-1,  launched  in  mid-1991,  will  have  a  part-time  scatterometer.  An 
identical  satellite,  ERS-2,  is  scheduled  for  1994.  The  U.S.  and  Japanese 
NSCAT/ADEOS  scatterometer  will  be  more  accurate,  and  it  will  sample 
continuously.  It  is  scheduled  for  1995.  The  ERS  satellites  will  both  carry 
altimeters,  which  could  yield  wind  data,  and  the  TOPEX/POSEIDON  al¬ 
timeter  is  scheduled  for  rnid-1992.  The  SMM/I  passive  microwave  radiometer 
carried  aboard  operational  Air  Force  DMSP  satellites  is  now  available  for  es¬ 
timating  wind  speeds,  and  this  series  is  likely  to  be  extended  in  the  future. 
Thus,  there  will  be  a  number  of  wind-sensing  satellites  available  in  the  next 
decade  which,  with  a  significant  effort,  could  be  combined  to  provide  oceanic 
wind  fields  of  greater  accuracy  than  available  now. 
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6  CONCLUSIONS  AND  RECOMMENDA¬ 
TIONS 

6.1  Conclusions 


Direct  flux  measurements  are  necessary  for  developing  correct  bulk  pa- 
rameterizations,  but  they  are  far  too  restricted  in  space  and  time  for  es¬ 
timating  global  average  CO2  fluxes.  Only  bulk  formul«is  are  practical  for 
estimating  space  and  time  averages  ais  <  >=<  E  •  ApC02  >.  In 

our  cursory  review  of  the  CO2  flux  problem,  we  find  important  uncertaini- 
ties  about:  1)  the  adequacy  of  ApC02  sampling  and  whether  the  data  are 
sufficient  for  computing  accurate  averages,  2)  the  dependence  of  E  on  wind 
speed  and  temperature,  3)  the  adequacy  of  climatological  winds  and  average 
ApC02  for  computing  average  fluxes,  and  4)  the  chemical  procedures  for 
calculating  ApC02.  In  view  of  these  uncertainities,  calculations  of  net  fluxes 
are  certain  to  disagree  so  wildly  that  we  conclude  it  is  not  known  how  to 
compute  the  present  air-to-ocean  flux  or  to  predict  future  changes.  Because 
predicting  future  CO2  concentrations  in  the  atmosphere  is  the  first  step  in 
assessing  global  warming  and  because  tramsport  of  CO2  to  the  ocean  is  com¬ 
parable  with  the  rate  of  accumulation  of  CO2,  we  perceive  an  urgent  need  to 
resolve  these  uncertaiinities  about  the  air-sea  CO2  flux. 

Both  ApC02  and  E  (which  depends  primarily  on  wind  speed)  vary  sub¬ 
stantially  in  time  and  space.  Although  extensive  measurements  of  ApC02 
have  been  made,  they  are  a  small  fraction  of  the  number  of  wind  speed  mea- 
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surements.  Furthermore,  only  a  small  fraction  of  the  collected  data  have 
been  published,  leading  to  uncertainities  about  the  quality  and  coverage  of 
ApC02  data  used  in  calculations.  Ov/ing  to  variations  of  CO2  concentration, 
temperature,  alkalinity,  and  total  dissolved  carbon  dioxide,  ApCOa  varies 
substantially  in  space  (see  Figures  3  and  4)  and  in  time.  Because  the  global 
average  of  ApCOa  is  only  a  few  ppm  (see  Figure  8  which  shows  zonal  av¬ 
erages)  and  much  less  than  its  variability,  accurate  space  and  time  averages 
are  essential. 

Few  repeated  observations  exist  for  determining  time  variability  in  PCO2, 
or  perhaps  more  importantly,  for  determining  the  factors  causing  the  vari¬ 
ability.  Wong  et  al.  [28]  and  collaborators  reportedly  have  analyzed  data 
from  Ocean  Weather  Station  P,  in  the  Gulf  of  Alaska,  and  found  that  ran¬ 
dom  variability  makes  it  impossible  to  detect  a  seasonal  cycle  or  determine 
an  accurate  mean  value.  If  this  is  generally  true,  a  time  average  can  cer¬ 
tainly  not  be  computed  from  two  or  three  samples.  Figures  3  and  4  show 
that  ApC02  is  subject  to  rapid,  front-like  variations  which  are  several  times 
larger  than  the  differences  in  the  zonal  averages  shown  in  Figure  11.  Conse¬ 
quently,  high  spatial  resolution  will  be  required  for  accurate  spatial  averages. 
At  the  very  least,  diligent  study  is  needed  to  define  the  sampling  needed  for 
realistic  space-time  averages. 

Figure  12  shows  the  ship  tracks  used  by  Tans  et  al.  (21).  Since  many 
of  the  data  were  obtained  from  commercial  ships  following  trade  routes,  cov¬ 
erage  is  good  in  northern  temperate  latitudes,  particularly  in  the  North  At¬ 
lantic.  Winter  observations  in  high  latitudes  are  rare,  even  in  the  North 
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Figure  12.  Ship  tracks  along  which  ApCO^  data  used  by  [Reference  21]  were  collected 

Compare  the  density  of  tracks  along  shipping  routes  in  the  North  and  South  Pacific 
with  the  Southern  Ocean  in  the  Pacific  and  Indian  sectors. 


Pacific,  but  especially  in  the  Southern  Ocean,  where  very  large  fluxes  are 
expected  in  winter.  Figure  3  from  Tans  et  al.  reveals  no  data  at  all  in  the 
Pacific  and  Indian  Ocean  sectors  of  the  Southern  Ocean  during  austral  win¬ 
ter.  With  data  gaps  of  this  magnitude,  it  is  hard  to  take  the  results  of  Tans 
et  al.  seriously. 

Although  the  Liss  and  Merlivat  [11]  parameterization  for  Vp  seems  to 
be  the  most  re^distic  at  present,  the  measurements  on  which  it  is  based  come 
mainly  from  laboratory  wind  tunnels  and  are  unlikely  to  represent  oceanic 
conditions  where  fluxes  are  expected  to  be  largest.  This  parameterization 
should  be  viewed  as  a  promising  beginning  rather  than  eis  a  definitive  result. 

Tans  et  al.  [21]  estimated  the  global  mean  flux,  <  Fcoj  >=<  E{U\q)  ■ 
ApC02  >  by  multiplying  averages  of  the  mean  wind  speeds  and  mean 
ApC02,  i.e.,  they  computed  <  Fcoj  >  =  ^(<  ^lo  >)•  <  ApC02  >.  The 

Vp  algorithm  of  Liss  and  Merlivat  [llj  shows  that  E  increases  nonlinearly  with 
wind  speed.  Consequently,  as  shown  by  Etcheto  and  Merlivat  [6]  using  cli¬ 
matological  winds  is  inadequate.  If  ApC02  were  constant  or  varied  indepen¬ 
dently  from  f/io,  then  it  would  follow  that  <  Fcoj=  F(<  U\q  >)•  <  ApC02  > 
and  <  E  >  could  be  computed  from  the  probability  density  of  Uio,  or  possi¬ 
bly  from  climatological  winds  by  applying  corrections  based  on  typical  fluc¬ 
tuations.  But  Tans  et  al.  did  not  account  for  fluctuations  and,  indeed,  in 
footnote  23  they  state  that  using  instantaneous  winds  with  Vp  from  Liss  and 
Merlivat  would  increase  absorption  by  25%  in  the  Northern  hemisphere  and 
50%  in  the  Southern  Ocean. 

While  independence  of  ApC02  and  Uiq  fluctuations  would  simplify  de- 


35 


termination  of  the  average  flux,  there  are  two  competing  reasons  why  there 
may  be  a  correlation  of  these  fluctuations.  First,  because  CO2  solubility  de¬ 
pends  on  temperature,  (pC02)oce  has  a  strong  temperature  dependence,  even 
if  the  concentration  [CO2]  remains  constant.  Strong  winds  deepen  the  mixed 
layer,  entraining  cold  water  into  it.  The  resulting  correlation  of  ApC02  and 
Uio  tends  to  increase  fluxes  into  the  ocean  over  what  would  be  found  without 
such  correlation.  At  the  same  time,  entrainment  typically  increases  [CO2]  in 
the  mixed  layer,  owing  to  higher  concentrations  of  CO2  in  the  deeper  water. 
Since  cooling  by  1°C  can  lower  (pC02)o<»  about  10  ppm,  it  is  quite  possible 
for  these  effects  to  change  the  sign  of  the  mean  flux  as  well  as  its  magnitude. 
For  example.  Figure  8  shows  ApC02  ^  0  over  much  of  the  North  Pacific,  so 
that  accounting  for  a  correlation  between  ApCO^  and  (/jo  nnay  significantly 
change  the  computed  flux.  Finally,  wind-driven  entrainment  also  brings  nu¬ 
trients  into  the  mixed  layer,  leading  to  the  possibility  that  biological  activity 
may  also  cause  additional  correlation  of  ApC02  and  wind  speed. 

Because  the  global  mean  ApC02  is  so  small,  it  is  critical  that  the  chem¬ 
ical  procedures  used  to  measure  it  be  accurate  to  a  very  few  ppm.  Although 
a  sample-to-sample  noise  several  times  larger  would  little  affect  the  averages, 
it  is  critical  that  there  not  be  a  mean  bias  error.  We  are  unaware  of  successful 
comparisons  of  different  methods  to  show  that  consistent  measurements  can 
be  made. 


36 


6.2  Recommendations 


In  view  of  the  wide  disparity  in  the  algorithms  and  in  the  data  used  for 
computing  CO2  fluxes,  little  can  be  expected  from  further  calculations  until 
better  algorithms  are  developed  and  the  data  improved.  Several  steps  are 
needed. 

1.  A  program  should  be  established  to  determine  Vp  at  sea,  as  is  done 
for  transfer  coefficients  for  momentum,  heat,  and  moisture.  This  will 
require  extensive  measurements  by  geochemists  and  turbulence  inves¬ 
tigators  working  at  the  same  place  and  time  from  platforms  and  ships. 
After  rationalizing  bulk  and  turbulent  fluxes  at  low  and  moderate  wind 
speeds,  the  measuremc  its  should  be  extended  to  the  high  winds  pro¬ 
ducing  the  hig^  g  s  transfer  rates  in  the  Southern  Ocean. 

2.  The  observational  work  to  determine  should  be  coupled  with  mea¬ 
surements  of  wave  breaking  and  subsurface  bubbles  and  with  theoreti¬ 
cal  and  laboratory  work  on  the  role  of  bubbles  in  gas  transfer. 

3.  Investigators  measuring  CO2  at  sea  should  be  funded  to  publish  their 
observations,  and  they  should  be  encouraged  to  adopt  common  meth¬ 
ods  of  measurement,  with  suitable  controls.  A  standard  is  needed  be¬ 
fore  beginning  systemmatic  sampling  of  ApC02. 

4.  Systenunatic  observations  should  be  established  to  obtain  adequate 
global  and  seasonal  coverage  of  pCO^-  At  present  observations  are 
constrained  to  piggyback  on  ship  tracks  established  for  other  purposes, 
leading  to  relatively  dense  sampling  of  the  North  Atlantic  and  sparse 
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sampling  of  the  Southern  Ocean.  Owing  to  the  fundamental  importance 
of  CO2  to  the  assessment  of  global  warming,  the  measurement  program 
should  seek  adequate  coverage  within  10  to  15  years. 

5.  Time  series  stations  should  be  established  to  determine  the  time  vari¬ 
ability  of  ApCOa  and  its  causes.  The  observations  should  include  mixed 
layer  parameters  to  test  whether  entrainment  causes  ApC02  to  vary 
and  whether  the  variation  can  be  parameterized  in  terms  of  wind  speed. 
Biological  parameters  should  also  be  monitored  to  assess  the  impor¬ 
tance  of  biological  productivity.  Similar  studies  may  be  needed  in  sev¬ 
eral  oceanographic  and  biological  regimes.  This  observational  program 
should  lead  a  vigorous  theoretical  program  to  interpret  the  observed 
variations  and  model  mean  fluxes  more  accurately  than  by  combining 
mean  ApC02  fields  and  mean  winds. 

6.  The  observational  base  of  ApC02  should  be  expanded  and  a  concerted 
effort  should  be  made  to  obtain  all  available  wind  observations  so  that 
the  global  average  flux  can  credible  estimated.  Analysis  of  time  series 
and  repeated  sections  should  be  used  to  estimate  the  sampling  errors 
associated  with  different  sampling  strategies  and  to  select  efficient  ones. 
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